Models of wealth distributions: a perspective 
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I. ABSTRACT 

A class of conserved models of wealth distributions are 
studied where wealth (or money) is assumed to be ex- 
changed between a pair of agents in a population just 
like the elastically colliding molecules of a gas exchang- 
ing energy. All sorts of distributions from exponential 
(Boltzmann-Gibbs) to something like Gamma distribu- 
tions and to that of Pareto's law (power law) are obtained 
out of such models with simple algorithmic exchange pro- 
cesses. Numerical investigations, analysis through tran- 
sition matrix and a mean field approach are employed to 
understand the generative mechanisms. A general sce- 
nario is examined wherefrom a power law and other dis- 
tributions can emerge. 



II. INTRODUCTION 

Wealth is accumulated in various forms and factors. 
The continual exchange of wealth (a value assigned) 
among the agents in an economy gives rise to interest- 
ing and many often universal statistical distributions of 
individual wealth. Here the word 'wealth' is used in a 
general sense for the purpose and the spirit of the re- 
view (inspite of separate meanings attached to the terms 
'money', 'wealth' and 'income'). Econophysics of wealth 
distributions [jj] is an emerging area where mainly the 
ideas and techniques of statistical physics are used in in- 
terpreting real economic data of wealth (available mostly 
in terms of income) of all kinds of people or other en- 
tities (e.g., companies) for that matter, pertaining to 
different societies and nations. Literature in this area 
is growing steadily (see an excellent website 01)- The 
prevalence of income data and apparent interconnections 
of many socio-economic problems with physics have in- 
spired a wide variety of statistical models, data analy- 
sis and other works in econophysics [3J , sociophysics and 
other emerging areas |^ in the last one decade or more 
(see an enjoyable article by Stauffer Q). 

Simple approach of agent based models have been able 
to bring out all kinds of wealth distributions that open 
up a whole new way of understanding and interpreting 
empirical data. One of the most important and contro- 
versial issues has been to understand the emergence of 
Pareto's law: 



where w > wq, wq being some value of wealth beyond 
which the power law is observed (usually towards the 
tail of the distributions). Pareto's law has been observed 
in income distributions among the people of almost all 
kinds of social systems across the world in a robust way. 
This phenomenon is now known for more than a century 
and has been discussed at a great length in innumer- 
able works in economics, econophysics, sociophysics and 
physics dealing with power law distributions. In many 
places, while mentioning Pareto's law, the power law is 
often written in the form: P{w) ex w^^^'^'^' , where f is 
referred to as 'Pareto index'. This index is usually found 
between 1 and 2 from empirical data fitting. Power laws 
in distributions appear in many other cases [a S 111 like 
that of computer file sizes, the growth of sizes of business 
firms and cities etc. Distributions are often referred to as 
'heavy tailed' or 'fat tailed' distributions ^]. Smaller the 
value of a, fatter the tail of the distribution as it may eas- 
ily be understood (the distribution is more spread out). 

Some early attempts 10] have been made to under- 
stand the income distributions which follow Pareto's law 
at the tail of the distributions. Some of them are stochas- 
tic logistic equations or some related generalized versions 
of that which have been able to generate power laws. 
However, the absence of direct interactions of one agent 
with any other often carries less significance in terms of 
interpreting real data. 

Some part of this review is centered around the concept 
of emergence of Pareto's law in the wealth distributions, 
especially in the context of the models that are discussed 
here. However, a word of caution is in order. In the cur- 
rent literature and as well as in the historical occurrences, 
the power law distribution has often been disputed with 
a closely related lognormal distribution |^. It is often 
not easy to distinguish between the two. Thus a brief 
discussion is made here on this issue. Let us consider the 
probability density function of a lognormal distribution: 
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If now the variance a^ in the lognormal distribution is 
large enough, the last term on the right hand side can 
be very small so that the distribution may appear linear 



on a log- log plot. Thus the cause of concern remains, 
particularly when one deals with real data. 

In the literature, sometimes one calculates a cumula- 
tive distribution function (to show the power law in a 
more convincing way) instead of plotting ordinary distri- 
bution from simple histogram (probability density func- 
tion). The cumulative probability distribution function 
P(> w) is such that the argument has a value greater 
than or equal to w: 
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P{> w)^ P{w')dw'. (4) 

J w 

If the distribution of data follows a power law P{w) — 
Cw~°', then 
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When the ordinary distribution (found from just his- 
togram and binning) is a power law, the cumulative dis- 
tribution thus also follows a power law with the exponent 
1 less: a — 1, which can be seen from a log- log plot of 
data. An extensive discussion on power laws and related 
matters can be found in Q • 

Besides power laws, a wide variety of wealth distri- 
butions from exponential to something like Gamma dis- 
tributions are all reported in recent literature in econo- 
physics. Exchange of wealth is considered to be a pri- 
mary mechanism behind all such distributions. In a 
class of wealth exchange models [ll|, [13, u^ Hjl that fol- 
low, the economic activities among agents have been as- 
sumed to be analogous to random elastic collisions among 
molecules as considered in kinetic gas theory in statisti- 
cal physics. Analogy is drawn between wealth (w) and 
Energy (E), where the average individual wealth (uJ) at 
equilibrium is equivalent to temperature (T). Wealth 
(w) is assumed to be exchanged between two randomly 
selected economic agents like the exchange of energy be- 
tween a pair of molecules in kinetic gas theory. The in- 
teraction is such that one agent wins and the other loses 
the same amount so that the sum of their wealth re- 
mains constant before and after an interaction (trading) : 
Wi{t+1)+Wj{t+1) = Wi{t)+Wj{t); each trading increases 
time t by one unit. Therefore, it is basically a process 
of zero sum exchange between a pair of agents; amount 
won by one agent is equal to the amount lost by another. 
This way wealth is assumed to be redistributed among 
a fixed number of agents (N) and the local conservation 
ensures the total wealth {W = J^'^i) o^ ^^^ the agents to 
remain conserved. 

Random exchange of wealth between a randomly se- 
lected pair of agents may be viewed as a gambling pro- 
cess (with zero sum exchange) which leads to Boltzmann- 
Gibbs type exponential distribution in individual wealth 
{P{w) ex exp(— w/uJ). However, a little variation in the 
mode of wealth exchange can lead to a distribution dis- 
tinctly different from ex pon ential. A number of agen t 
based conserved models 



invoked in recent times, are essentially variants of a gam- 
bling process. A wide variety of distributions evolve out 
of these models. There has been a renewed interest in 
such two-agent exchange models in the present scenario 
while dealing with various problems in social systems 
involving complex interactions. A good insight can be 
drawn by looking at the 2x2 transition matrices associ- 
ated with the process of wealth exchange 22] . 

In this review, the aim would be to arrive at some un- 
derstanding of how wealth exchange processes in a sim- 
ple working way may lead to a variety of distributions 
within the framework of the conserved models A fixed 
number of TV agents in a system are allowed to interact 
(trade) stochastically and thus wealth is exchanged be- 
tween them. The basic steps of such a wealth exchange 
model are as follows: 



w^{t + l) = w^{t)-\-Aw, 



Wj{t + 1) = Wj{t) — Aw, 



(6) 
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where Wi (t) and Wj (t) are wealths of i-th and j-th agents 
at time t and Wi {t -\- 1) and Wj (t-\-l) are that at the next 
time step {t-\-l). The amount Aw (to be won or to be lost 
by an agent) is determined by the nature of interaction. 
If the agents are allowed to interact for a long enough 
time, a steady state equilibrium distribution for individ- 
ual wealth is achieved. The equilibrium distribution does 
not depend on the initial configuration (initial distribu- 
tion of wealth among the agents). A single interaction 
between a randomly chosen pair of agents is referred here 
as one 'time step'. In some simulations, N such interac- 
tions are considered as one time step. This, however, 
does not matter as long as the system is evolved through 
enough time steps to come to a steady state and then 
data is collected for making equilibrium probability dis- 
tributions. For all the numerical results presented here, 
data have been produced following the available mod- 
els, conjectures and conclusions. Systems of TV = 1000 
agents have been considered in each case. In each numer- 
ical investigation, the system is allowed to equilibrate for 
a sufficient time that ranges between 10^ to 10* time 
steps. Configuration averaging has been done over 10'^ 
to 10^ initial configurations in most cases. The average 
wealth (averaged over the agents) is kept fixed at w = 1 
(by taking total wealth, W — N) for all the cases. The 
wealth distributions, that are dealt here in this review, 
are ordinary distributions (probability density function) 
and not the cumulative ones. 



III. PURE GAMBLING 

In a pure gambling process (usual kinetic gas theory) , 
entire sum of wealths of two agents is up for gambling. 
Some random fraction of this sum is shared by one agent 
and the rest goes to the other. The randomness or 



stochasticity is introduced into the model through a pa- 
rameter e which is a random number drawn from a uni- 
form distribution in [0, 1]. (Note that e is independent 
of a pair of agents i.e., a pair of agents is not hkely to 
share the same fraction of aggregate weahh the same way 
when they interact repeatedly). The interaction can be 
seen through: 



w.,it + 1) = e[w,{t) + Wj{t)], 



Wj{t + 1) = {1 - e)[w^{t) + Wj{t)], 



(7) 



where the pair of agents (indicated by i and j) are cho- 
sen randomly. The amount of wealth that is exchanged 
is Aw = e[wi{t) + Wj{t)] — Wi{t). The individual wealth 
distribution {P{w) vs. w) at equilibrium emerges out to 
be Boltzmann-Gibbs distribution like exponential. Expo- 
nential distribution of personal income has in fact been 
shown to appear in real data jla Il4| . In the kinetic 
theory model, the exponential distribution is found by 
standard formulation of master equation or by entropy 
maximization method, the latter has been discussed later 
in brief in section rVIIII A normalized exponential distri- 
bution obtained numerically out of this pure gambling 
process is shown in Fig. ^ in semi logarithmic plot. The 
high end of the distribution appears noisy due to sam- 
pling of data. The successive bins on the right hand side 
of the graph contain less and less number of samples in 
them so the fractional counts in them are seen to fluctu- 
ate more (finite size effect). One way to get rid of this 
sampling error in a great extent is by way of taking log- 
arithmic binning Q. Here it is not important to do so 
as the idea is to show the nature of the curve only. (In 
the case of power law distribution, an even better way 
to demonstrate and extract the power law exponent is to 
plot the cumulative distribution as discussed already.) 
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FIG. 1: Distribution of wealth for the case of Pure Gambling: 
the linearity in the semi-log plot indicates exponential distri- 
bution. 

If one intends to take time average of wealth of a sin- 
gle agent over a sufficient time, it comes out to be equal 



for all the agents. Therefore, the distribution of individ- 
ual time averaged wealth turns out to be a delta function 
which is checked from numerical data. This is because 
the fluctuation of wealth of any agent over time is sta- 
tistically no different from that of any other. The same 
is true in case of the distribution of wealth of a single 
agent over time. However, when the average of wealth of 
any agent is calculated over a short time period, the delta 
function broadens and the right end part of which decays 
exponentially. The distribution of individual wealth at a 
certain time turns out to be purely exponential as men- 
tioned earlier. This may be thought of as a 'snap shot' 
distribution. 



IV. UNIFORM SAVING PROPENSITY 

Instead of random sharing of their aggregate wealth 
during each interaction, if the agents decide to save (keep 
aside) a uniform (and fixed) fraction (A) of their current 
individual wealth, then the wealth exchange equations 
look like the following: 

w^it + 1) ^ Xw,{t) + e{l - X) [w, {t) + wj {t)] , (8) 

Wj{t + 1) = Xwjit) + (1 - e)(l - X)[wi{t) + Wj{t)], 

where the amount of wealth that is exchanged is Aw = 
(e— l){l — X)[wi{t) +Wj{t)\. The concept of saving as in- 
troduced by Chakrabarti and group |l2| | in an otherwise 
gambling kind of interactions brings out distinctly differ- 
ent distributions. A number of numerical works followed 
123, U% I23 in order to understand the emerging distri- 
butions to some extent. Saving induces accumulation of 
wealth. Therefore, it is expected that the probability of 
finding agents with zero wealth may be zero unlike in 
the previous case of pure gambling where due to the un- 
guarded nature of exchange many agents are likely to go 
nearly bankrupt! (It is to be noted that for an exponen- 
tial distribution, the peak is at zero.) In this case the 
most probable value of the distribution (peak) is some- 
where else than at zero (the distribution is right skewed) . 
The right end, however, decays exponentially for large 
values of w. It has been claimed through heuristic argu- 
ments (based on numerical results) that the distribution 
is a close approximate form of the Gamma distribution 
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where the Gamma function T{n) and the index n are un- 
derstood to be related to the saving propensity parameter 
A through the following relation: 



1 
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(10) 



The emergence of probable Gamma distribution is also 
subsequently supported through numerical results in |24| . 



However, it has later been shown in [2g| by considering 
moments' equation that moments up to third order agree 
with that obtained from the above form of distribution 
subject to the condition stated in eqn. (|10|l . Discrepan- 
cies start showing only from 4th order onwards. There- 
fore, the actual form of distribution still remains an open 
question. 

In Fig. [21 two distributions are shown for two different 
values of saving propensity factor: A = 0.4 and A = 0.8. 
Smaller the value of A, lesser the amount one is able to 
save. This in turn means more wealth is available in the 
market for gambling. In the limit of zero saving (A = 0) 
the model reduces to that of pure gambling. In the oppo- 
site extent of large saving, only a small amount of wealth 
is up for gambling. Then the exchange of wealth will not 
be able to drastically change the amount of individual 
wealth. This means the width of distribution of individ- 
ual wealth will be narrow. In the limit of A = 1, all the 
agents save all of their wealth and thus the distribution 
never evolves. The concept of 'saving' here is of course 
a little different from that in real life where people do 
save some amount to be kept in a bank or so and the in- 
vestment (or business or gambling) is done generally not 
with the entire amount (or a fraction) of wealth that one 
holds at a time. 



Wj{t + l)^ Xwj{t). 

The above equations indicate that the randomly selected 
agent j keeps (saves) an amount Xwj (t) which is propor- 
tional to the wealth he currently has and transfers the 
rest to the other agent i. This is indeed a stochastic pro- 
cess and is able to produce Gamma type distributions in 
wealth as observed. However, a distribution with ran- 
dom e and that with a fixed e are different. Numerically, 
it has been observed, the distribution with A = 0.8 and 
with e = 1 is very close to that with A = 0.5 and with 
random e. In Fig. |21 the distribution with fixed A — 0.8 
and fixed e = 1 is plotted along with other two distribu- 
tions with random e. It should also be noted that while 
with fixed e, one does not get Gamma type distributions 
for all values of A; especially for low values of A the dis- 
tributions become close to exponential as observed. This 
is not clearly understood though. 

It has recently been brought to notice in [23 that a 
very similar kind of agent based model was proposed by 
Angle [23 (see other references cited in |23) in sociolog- 
ical journals quite some years ago. The pair of equations 
in Angle's model are as follows: 

w^it + 1) = w,{t) + DtLuWjit) - (1 - Dt)ujw^{t), (12) 
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FIG. 2: Wealth distribution for the model of uniform and 
fixed saving propensity. Two distributions are shown with 
A — 0.4 and A — 0.8 where the stochasticity parameter e is 
drawn randomly and uniformly in [0, 1]. Another distribution 
is plotted with A = 0.8 but with fixed value of the stochasticity 
parameter, e — 1. 



Stochastic evolution of individual wealth is also exam- 
ined without the inclusion of stochastic parameter e. The 
stochasticity seems to be automatically introduced any- 
how through the random selection of a pair of agents (and 
the random choice of the winner or loser as well) at each 
time. Therefore, it is interesting to see how the distribu- 
tions evolve with a fixed value of e. As an example, the 
equations in JSj) reduce to the following with e = 1 : 



W^{t + 1) = W^{t) + (1 - X)Wj{t), 



(11) 



Wj{t + 1) = Wj{t) + (1 - Dt)uJWi{t) - DtUJWj{t), 

where w is a fixed fraction and the winner is decided 
through a random toss Dt which takes a value either or 
1. Now, the above can be seen as the more formal way of 
writing the pair of equations (|11|) which can be arrived 
at by choosing Dt = 1 and identifying w = (1 — A). 

It can in general be said, within the framework of this 
kind of (conserved) models, different ways of incorpo- 
rating wealth exchange processes may lead to drastically 
different distributions. If the gamble is played in a biased 
way, then this may lead to a distinctly different situation 
than the case when it is played in a normal unbiased 
manner. Since in this class of models negative wealth or 
debt is not allowed, it is desirable that in each wealth ex- 
change, the maximum that any agent may invest is the 
amount that he has at that time. Suppose, the norm 
is set for an 'equal amount invest' where the amount to 
be deposited by an agent for gambling is decided by the 
amount the poorer agent can afford and consequently 
the same amount is agreed to be deposited by the richer 
agent. Let us suppose Wi > Wj. Now, the poorer agent 
(j) may invest a certain fraction of his wealth, an amount 
Xwj and the rest {1—X)wj is saved by him. Then the total 
amount 2Xwj is up for gambling and as usual a fraction of 
this, 2eXwj may be shared by the richer agent i where the 
rest (1 — e)Xwj goes to the poorer agent j. This may ap- 
pear fine, however, this leads to 'rich gets richer and poor 
gets poorer' way. The richer agent draws more and more 
wealth in his favour in the successive encounters and the 
poorer agents are only able to save less and less and fi- 
nally there is a condensation of wealth at the hands of the 



richest person. This is more apparent when one considers 
an agent with A = 1 where it can be easily checked that 
the richer agent automatically saves an amount equal to 
the difference of their wealth (wt — Wj) and the poorer 
agent ends up saving zero amount. Eventually, poorer 
agents get extinct. This is 'minimum exchange model' 



V. DISTRIBUTED SAVING PROPENSITY 

The distributions emerge out to be dramatically differ- 
ent when the saving propensity factor (A) is drawn from a 
uniform and random distribution in [0,1] as introduced in 
a model proposed by Chatterjee, Chakrabarti and Manna 
[l5|. Randomness in A is assumed to be quenched {i.e., 
remains unaltered in time). Agents are indeed heteroge- 
neous. They are likely to have different (characteristic) 
saving propensities. The pair of wealth exchange equa- 
tions are now written as: 

w,it+l) = X,w,{t)+e[il-X,)w,it) + {l~Xj)wjit))], (13) 

Wjit + 1) = X,w,it) + {l-e)[{l-X,)w^it) + il-\,)w,m. 

A power law with exponent a = 2 (Pareto index v = 1) 
is observed at the right end of the wealth distribution for 
several decades. Such a distribution is plotted in Fig. |2| 
where a straight line is drawn in the log-log plot with 
slope = -2 to illustrate the power law and the exponent. 
Extensive numerical results with different distributions in 
the saving propensity parameter A are reported in [l8l| . 
Power law (with exponent a = 2) is found to be ro- 
bust. The value of Pareto index obtained here {ly — 1), 
however, differs from what is generally extracted (1.5 or 
above) from most of the empirical data of income dis- 
tributions (see discussions and analysis on real data by 
various authors in [ij). The present model is not able 
to resolve this discrepancy and it is not expected at the 
outset either, introducing random waiting time in the 
interactions of agents in order to have a justification for 
a larger value of the exponent v {2 ft] . 

The distributed saving gives rise to an additional inter- 
esting feature when a special case is considered where the 
saving parameter A is assumed to have only two fixed val- 
ues, Ai and A2 (preferably widely separated) . A bimodal 
distribution in individual wealth results in |22l | . This can 
be seen from the Fig. ^ The system evolves towards a 
robust and distinct two-peak distribution as the differ- 
ence in Ai and A2 is increased systematically. Later it is 
seen that one still gets a two-peak distribution even when 
Ai and A2 are drawn from narrow distributions centered 
around two widely separated values (one large and one 
small). Two economic classes seem to persist until the 
distributions in Ai and A2 have got sufficient widths. A 
population can be imagined to have two distinctly differ- 
ent kinds of people: some of them tend to save a very 
large fraction (fixed) of their wealth and the others tend 
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FIG. 3: Wealth distribution for the model of random saving 
propensity plotted in log-log scale. A straight line with slope 
= -2 is drawn to demonstrate that the power law exponent is 
a = 2. 
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FIG. 4: Bimodal distribution of wealth (w) with fixed values 
of saving propensities, Ai=0.2 and A2=0.8. Emergence of two 
economic classes are apparent. 



to save a relatively small fraction (fixed) of their wealth. 
Bimodal distributions (and a polymodal distribution, in 
general) are, in fact, reported with real data for the in- 
come distributions in Argentina 30] . The distributions 
were derived at a time of political crisis and thus they 
may not be regarded as truly equilibrium distributions 
though. However, it remains an interesting possibility 
out of a simple model of wealth exchange. 



A. Po\ver law from mean field analysis 

One can have an estimate of ensemble averaged value 
of wealth JU using one of the above equations itT^ in 
section^ Emergence of a power law in the wealth distri- 
bution can be established through a simple consideration 
as follows. Taking ensemble average of all the terms on 



both sides of the first eqn. (|13|l . one may write: 



1 



N 






The last term on the right hand side is replaced by the 
average over agents where it is assumed that any agent 
(here the i-th agent), on an average, interacts with all 
other agents of the system, allowing sufficient time to 
interact. This is basically a mean field approach. If e 
is assumed to be distributed randomly and uniformly 
between and 1 then (e) — i. Wealth of each indi- 
vidual keeps on changing due to interactions (or wealth 
exchange processes that take place in a society) . No mat- 
ter what the personal wealth one begins with, the time 
evolution of wealth of an individual agent at the steady 
state is independent of that initial value. This means 
the distribution of wealth of a single agent over time is 
stationary. Therefore, the time averaged value of wealth 
of any agent remains unchanged whatever the amount 
of wealth one starts with. In course of time, an agent 
interacts with all other agents (presumably repeatedly) 
given sufficient time. One can thus think of a number 
of ensembles (configurations) and focus attention on a 
particular tagged agent who eventually tends to interact 
with all other agents in different ensembles. Thus the 
time averaged value of wealth is equal to the ensemble 
averaged value in the steady state. 
Now if one writes 



1 ^ 

((1-A).«)^(-5](1-A. 



i=i 



the above equation H14|l reduces to: 



il-X,){w,) = {{l-X)w)), 



(15) 



(16) 



The right hand side of the above equation is independent 
of any agent-index and the left hand side is referred to 
any arbitrarily chosen agent i. Thus, it can be argued 
that the above relation can be true for any agent (for 
any value of the index i) and so it can be equated to a 
constant. Let us now recognize C = ((1 — X)w), a con- 
stant which is found by averaging over all the agents in 
the system and which is further averaged over ensembles. 
Therefore, one arrives at a unique relation for this model: 



C 



(1-A)' 



(17) 



where one can get rid of the index i and may write 
(wi) = w for brevity. The above relation is also veri- 
fied numerically which is obtained by many authors in 
their numerical simulations and scaling of data [H, |2J| . 

2 

One can now derive dw = ^dX from the above relation 
(|17|l . An agent with a (characteristic) saving propensity 



factor (A) ends up with wealth (w) such that one can in 
general relate the distributions of the two: 



P{w)dw = g{X)dX. 



(18) 



If now the distribution in A is considered to be uniform 
then 5(A) — constant. Therefore, the distribution in w is 
bound to be of the form: 



p{w) (X 



1 



w^ 



(19) 



This may be regarded as Pareto's law with index a — 
2 which is already numerically verified for this present 
model. The same result is also obtained recently in |23| 
where the treatment is argued to be exact. 



B. Power law from reduced situation 

From numerical investigations, it seems that the 
stochasticity parameter e is irrelevant as long as the sav- 
ing propensity parameter A is made random. It has been 
tested that the model is still able to produce power law 
(with the same exponent, a — i) for any fixed value of e. 
As an example, the case for e = 1 is considered. The pair 
of wealth exchange equations (refakg:eqn:ransave) now 
reduce to the following: 

w,{t + l) = w^{t) + (1 - Xj)wj{t) ^ w^{t) + iqjWj{t), (20) 

Wj{t + 1) = w,{t) - (1 - Xj)wj{t) - (1 - fij)w,{t). 

The exchange amount. Aw = (1 — Xj)wj{t) — r]jWj{t) is 
now regulated by the parameter rj — (1 — A) only. If A is 
drawn from a uniform and random distribution in [0, 1], 
then r] is also uniformly and randomly distributed in [0, 
1]. To achieve a power law in the wealth distribution it 
seems essential that randomness in rj has to be quenched. 
For 'annealed' type disorder (i.e., when the distribution 
in r] varies with time) the power law gets washed away 
(which is observed through numerical simulations). It 
has also been observed that power law can be obtained 
when rj is uniformly distributed between and some value 
less than or equal to 1. As an example, rj is taken in the 
range between and 0.5, a power law is obtained with 
the exponent around a = 2. However, when 77 is taken in 
the range 0.5 < 77 < 1, the distribution clearly deviates 
from power law which is evident from the log-log plot in 
Fig. |5l Thus there seems to be a crossover from power 
law to some distribution with exponentially decaying tail 
as one tunes the range in the quenched parameter rj. 

At this point, two important criteria may be identified 
for achieving power law within this reduced situation: 

• The disorder in the controlling parameter 77 has to 
be quenched (fixed set of Ty's for a configuration of 
agents) , 

• It is required that 77, when drawn from a uniform 
distribution, the lower bound of that should be 0. 
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FIG. 5: Wealth distributions (plotted in the log-log scale) 
for two cases of the 'reduced situation': (i) < r; < 0.5 
and (ii) 0.5 < r] < 1 plotted in log-log scale. In one case, 
the distribution follows a power law (with exponent around 
a = 2) and in the other case, it is seen to be clearly deviating 
from a power law. 



The above criteria may appear ad hoc, nevertheless have 
been checked by extensive numerical investigations. It 
is further checked that the power law exponent does not 
depend on the width of the distribution in r] (as long as 
it is between and something less than 1). This claim 
is substantiated by taking various ranges of r] in which 
it is uniformly distributed. Systematic investigations are 
made for the cases where rj is drawn in [0, 0.2], [0, 0.4], 
... ,[0, 1]. Power laws result in in all the cases with the 
exponent around a = 2. 



VI. UNDERSTANDING THROUGH 
TRANSITION MATRIX 

The evolution of wealth in the kind of two-agent wealth 
exchange process can be described through the following 
2x2 transition matrix (T) 22]: 



= T 



W-i 



where it is written, w'^ = Wi(t + I) and Wi = Wi(t) and 
so on. The transition matrix (T) corresponding to pure 
gambling process (in section HTH) can be written as: 



T 



e e 

1-e 1-e 



In this case the above matrix is singular (determinant, 
\T\ = 0) which means the inverse of this matrix does not 
exit. This in turn indicates that an evolution through 
such transition matrices is bound to be irreversible. This 
property is connected to the emergence of exponential 
(Boltzmann-Gibbs) wealth distribution. The same may 



be perceived in a different way too. When a product of 
such matrices (for successive interactions) are taken, the 
left most matrix (of the product) itself returns: 



1-e 1-e 



1-ei 1-ei 



e e 

1-e 1-e 



The above signifies the fact that during the repeated 

interactions of the same two agents (via this kind 

of transition matrices), the last of the interactions is 

what matters (the last matrix of the product survives) 
[j.(n)_j.(n-i) ^ ^ ^ y(2) y(i) ^ j,(„)]^ rpj^jg .^Qgg Qf memory' 

(random history of collisions in case of molecules) may 
be attributed here to the path to irreversibility in time. 
The singularity can be avoided if one considers the 
following general form: 



Ti 



ei £2 

1 - ei 1 - e2 



where ei and 62 are two different random numbers drawn 
uniformly from [0, 1] (This ensures the transition matrix 
to be nonsingular.). The significance of this general form 
can be seen through the wealth exchange equations in the 
following way: the ei fraction of wealth of the 1st agent 
(i) added with £2 fraction of wealth of the 2nd agent (j) 
is retained by the 1st agent after the trade. The rest of 
their total wealth is shared by the 2nd agent. This may 
happen in a number of ways which can be related to the 
detail considerations of a model. The general matrix Ti 
is nonsingular as long as ei 7^ €2 and then the two-agent 
interaction process remains reversible in time. Therefore, 
it is expected to have a steady state equilibrium distribu- 
tion of wealth which may deviate from exponential distri- 
bution (as in the case with pure gambling model). When 
one considers ei — £2, one again gets back the pure expo- 
nential distribution. A trivial case is obtained for ei = 1 
and e2 = 0. The transition matrix then reduces to the 
■ 1 



identity matrix / = 



1 



which trivially corresponds 



for such conserved models is bound to be 



to no interaction and no evolution. 

It may be emphasized that any transition matrix 
ill ^12 

^21 ^22 

of the form such that the sum of two elements of either 
of the columns has to be unity by design: tn + t2i = 1, 
^12 + ^22 = 1- It is important to note that whatever ex- 
tra parameter, no matter, one incorporates within the 
framework of the conserved model, the transition matrix 
has to retain this property. 

In Fig. 13 three distributions (with ei ^ £2) are plot- 
ted where ei and €2 are drawn randomly from uniform 
distributions with different ranges. It is demonstrated 
that qualitatively different distributions are possible as 
the parameter ranges are tuned appropriately. 

Now let us compare the above situation with the model 
of equal saving propensity as discussed in section IIVI 



8 







0.5 1 1.5 2 2.5 3 3.5 4 
w (wealth) 

FIG. 6: Three normalized wealth distributions are shown cor- 
responding to the general matrix T2 (in general with ei 7^ £2) 
as discussed in the text. Curves are marked by numbers (1,2 
and 3) and the ranges of ei and £2 are indicated within which 
they are drawn uniformly and randomly. 



With the incorporation of saving propensity factor A, the 
transition matrix now looks like: 



A + e(l-A) e(l-A) 

(l-e)(l-A) A+(l-e)(l-A) 



The matrix elements can now be rescaled by assuming 
ei = A + e(l — A) and £2 = e(l — A) in the above matrix. 
Therefore, the above transition matrix reduces to 



T2 = 



£1 £2 

1 - ei 1-62 



Thus the matrix T2 is of the same form as Ti. The dis- 
tributions due to above two matrices of the same general 
form can now be compared if one can correctly identify 
the ranges of the rescaled elements. In the model of 
uniform saving: A < ei < 1 and < £2 < (1 — A) as 
the stochasticity parameter e is is drawn from a uniform 
and random distribution in [0, 1]. As long as ei and £2 
are different, the determinant of the matrix is nonzero 
(|2^2| = ei — £2 = A). Therefore, the incorporation of the 
saving propensity factor A brings two effects: 

• The transition matrix becomes nonsingular, 

• The matrix elements tn (— ei) and ti2 {— £2) are 
now drawn from truncated domains (somewhere in 
[0, 1]). 

Hence it is clear from the above discussion that the 
wealth distribution with uniform saving is likely to be 
qualitatively no different from what can be achieved with 
general transition matrices having different elements, 
ei 7^ £2- The distributions obtained with different A may 
correspond to that with appropriately chosen ei and £2 
in Ti. 



In the next stage, when the saving propensity factor 
A is distributed as in section ]V\ the transition matrix 
between any two agents having different A's (say, Ai and 
A2) now looks like: 



Ai+£(1-Ai) £(1-A2) 

(1-£)(1-Ai) A2 + (l-£)(l-A2) 



Again the elements of the above matrix can be rescaled 
by putting e'^ = Ai + £(1 — Ai) and £3 = £(1 — A2). Hence 
the transition matrix can again be reduced to the same 
form as that of Ti or T2: 



T. 



1 - e~ 



l-e'2 



The determinant here is jTaj = e[ — ^'2 = Ai(l — £) + £A2. 
Here also the determinant is ensured to be nonzero as all 
the parameters £, Ai and A2 are drawn from the same 
positive domain: [0, 1]. This means that each transition 
matrix for two-agent wealth exchange remains nonsingu- 
lar which ensures the interaction process to be reversible 
in time. Therefore, it is expected that qualitatively dif- 
ferent distributions are possible when one appropriately 
tunes the two independent elements in the general form 
of transition matrix (Ti or T2 or T3). However, the emer- 
gence of power law tail {Pareto's law) in the distribution 
can not be explained by this. Later it is examined that 
to obtain a power law in the framework of present mod- 
els, it is essential that the distribution in A has to be 
quenched (frozen in time) which means the matrix ele- 
ments in the general form of any transition matrix have 
to be quenched. In the section IV Bl it has been shown 
that the model of distributed saving (section^) is equiv- 
alent to a reduced situation where one needs only one 
variable 77. The this case the corresponding transition 
matrix looks even simpler: 



Ta 



1 j; 
1-7/ 



where a nonzero determinant (IT4I = 1 — 77 =/= 0) is ensured 
among other things. 



A. Distributions from generic situation 

From all the previous discussions, it is clear that the 
the transition matrix (for zero sum wealth exchange) is 
bound to be of the following general form: 



ei £2 

1 - £1 1 - £2 



The matrix elements, £1 and £2 can be appropriately as- 
sociated with the relevant parameters in a model. A 



generic situation arrives where one can generate all sorts 
of distributions by controlling ei and £2- 

As long as ei 7^ £2, the matrix remains nonsingular 
and one achieves Gamma type distributions. In a spe- 
cial case, when ei = £2, the transition matrix becomes 
singular and Boltzmann-Gibbs type exponential distri- 
bution results in. It has been numerically checked that a 
power law with exponent a = 2 is obtained with the gen- 
eral matrix when the elements £1 and £2 are of the same 
set of quenched random numbers drawn uniformly in [0, 
1]. The matrix corresponding to the reduced situation in 
the section IVBI as discussed, is just a special case with 
€1 = 1 and £2 = r], drawn from a uniform and (quenched) 
random distribution. Incorporation of any parameter in 
an actual model (saving propensity, for example) results 
in the adjustment or truncation of the full domain [0, 1] 
from which the element £1 or £2 is drawn. Incorporating 
distributed A's in section [3 is equivalent to considering 
the following domains: Ai < £1 < 1 and < £2 < (1 — A2). 

A more general situation arrives when the matrix ele- 
ments £1 and £2 are of two sets of random numbers drawn 
separately (one may identify them as ei and £3 to dis- 
tinguish) from two uniform and random distributions in 
the domain: [0, 1]. In this case a power law is obtained 
with the exponent a = 3 which is, however, distinctly dif- 
ferent from that obtained in 'distributed saving model' 
in section^ To test the robustness of the power law, the 
distributions in the matrix elements are taken in the fol- 
lowing truncated ranges: 0.5 < £1 < 1 and < £2 < 0.5 
(widths are narrowed down). A power law is still ob- 
tained with the same exponent {a close to 3). These 
results are plotted in Fig. [7| 
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FIG. 7: Distribution of individual wealth (w) for the most 
general case with random and quenched ei and £2: The 
elements are drawn from two separate distributions where 
< ei < 1 and < £2 < 1 in one case and in the other 
case, they are chosen from the ranges, 0.5 < ei < 1 and 
< £2 < 0.5. Both show power laws with the same expo- 
nent around 3.0 (the two distributions almost superpose). A 
straight line (with slope -3.0) is drawn to demonstrate the 
power law in the log-log scale. 

It is possible to achieve distributions other than power 



laws as one draws the matrix elements, £1 and £2 from dif- 
ferent domains within the range between and 1. There 
is indeed a crossover from power law to Gamma like dis- 
tributions as one tunes the elements. It appears from ex- 
tensive numerical simulations that power law disappears 
when both the parameters are drawn from some ranges 
that do not include the lower limit 0. For example, when 
it is considered, 0.8 < £1 < 1.0 and 0.2 < £2 < 0.4, 
the wealth distribution does not follow a power law. In 
contrast, when £1 and £2 are drawn from the ranges, 
0.8 < £1 < 1.0 and < £2 < 0.1, the power law dis- 
tribution is back again. 

It now appears that to achieve a power law in such a 
generic situation, the following criteria are to be fulfilled: 

• It is essential to have the randomness or disorder 
in the elements £1 and £2 to be quenched, 

• In the most general case, £1 should be drawn from 
a uniform distribution whose upper bound has to 
be 1 and for £2 the lower bound has to be 0. Then a 
power law with higher exponent a = 3 is achieved. 
To have a power law with exponent a = 2, the ma- 
trix elements are to be drawn from the same dis- 
tribution. (These choices automatically make the 
transition matrices to be nonsingular.) 

The above points are not supported analytically at this 
stage. However, the observation seems to bear important 
implications in terms of generation of power law distri- 
butions. 

When the disorder or randomness in the elements £1 
and £2 change with time {i.e., not quenched) unlike the 
situation just discussed above, the problem is perhaps 
similar to the mass diffusion and aggregation model by 
Majumdar, Krishnamurthy and Barma [SJI- The mass 
model is defined on a one dimensional lattice with pe- 
riodic boundary condition. A fraction of mass from a 
randomly chosen site is assumed to be continually trans- 
ported to any of its neighbouring sites at random. The 
total mass between the two sites then is unchanged (one 
site gains mass and the other loses the same amount) and 
thus the total mass of the system remains conserved. The 
mass of each site evolves as 



m,{t+ 1) = (1 - r]i)mi{t) + rijmj{t). 



(21) 



Here it is assumed that r]i fraction of mass mi is dissoci- 
ated from that site i and joins either of its neighbouring 
sites J = I ± 1. Thus (1 — 77^) fraction of mass mt re- 
mains at that site whereas a fraction rjj of mass mj from 
the neighbouring site joins the mass at site i. Now if we 
identify £1 = (1 — 77^) and £2 = rjj then this model is just 
the same as described by the general transition matrix 
as discussed so far. If jy^'s are drawn from a random and 
uniform distribution in [0, 1] then a mean field calcula- 
tion (which turns out to be exact in the thermodynamic 
limit), as shown in i3^, brings out the stationary mass 
distribution P{m) to be a Gamma distribution: 



4?Tl 

P[ra) — z^^ 
m 



2m/m 



(22) 
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where rn is the average mass of the system. It has been 
numericaUy checked that there seems to be no apprecia- 
ble change in the distribution even when the lattice is not 
considered. Lattice seems to play no significant role in 
the case of kinetic theory like wealth distribution mod- 
els as well. Incidentally, this distribution with m = 1 is 
exactly the same as the Gamma distribution [eqn. (jS))], 
mentioned in section IIVI when one considers n = 2 in 
that. The index n equals to 2 when one puts A = ^ in 
the relation ((Tn|l . 

In the general situation (ei =/= £2), when both the pa- 
rameters are drawn from a random and uniform distribu- 
tion in [0, 1], the emerging distribution very nearly fol- 
lows the above expression (|22() . Only when the random- 
ness in them is quenched (fixed in time), there is a pos- 
sibility of getting a power law as it is already mentioned. 
The Gamma distribution [eqn. (|22|l ] and the numerically 
obtained distributions for different cases (as discussed in 
the text) are plotted in Fig. |S1 in support of the above 
discussions. 
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FIG. 8: Normalized probability distribution functions ob- 
tained for three different cases: (i) Wealth distribution with 
random and uniform ei and €2 in [0, 1], (ii) Wealth distri- 
bution with uniform and fixed saving propensity, A = |, (iii) 
Mass distribution for the model [3^ in one dimensional lattice 
(as discussed in text). The theoretical Gamma distribution 
[the eqn. 1221 1 is also plotted (line draw) to have a comparison. 



VII. ROLE OF SELECTIVE INTERACTION 

So far the models of wealth exchange processes have 
been considered where a pair of agents is selected ran- 
domly. However, interactions or trade among agents in 
a society are often guided by personal choice or some 
social norms or some other reasons. Agents may like to 
interact selectively and it would be interesting to see how 
the Individual wealth distribution is influenced by selec- 
tion [13 . The concept of selective interaction is already 
there when one considers the formation of a family. The 
members of a same family are unlikely to trade (or inter- 



act) among each other. It may be worth to examine the 
role played by the concept of 'family' in wealth distribu- 
tions of families: 'family wealth distribution' for brevity. 
A family in a society usually consists of more than one 
agent. In computer simulation, the agents belonging to a 
same family are coloured to keep track of. To find wealth 
distributions of families, the contributions of the same 
family members are added up. In Fig. |^ family wealth 
distributions are plotted for three cases: (i) families con- 
sist of 2 members each, (ii) families consist of 4 members 
each and (iii) families of mixed sizes between 1 and 4. 
The distributions are clearly not pure exponential, but 
modified exponential distributions (Gamma type distri- 
butions) with different peaks and different widths. This 
is quite expected as the probability of zero income of a 
family is zero. Modified exponential distribution for fam- 
ily wealth is also supported by fitting real data [l4J . 
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FIG. 9: Family wealth distributions: two curves are for fam- 
ilies consisting of all equal sizes of 2 and 4. One curve is for 
a system of families consisting of various sizes between 1 and 
4. The distributions are not normalized. 

Some special way of incorporating selective interaction 
is seen to have a drastic effect in the individual wealth dis- 
tribution. To implement the idea of 'selection', a 'class' 
of an agent is defined through an index e. The class may 
be understood in terms of some sort of efficiency of accu- 
mulating wealth or some other closely related property. 
Therefore, e's are assumed to be quenched. It is assumed 
that during the interactions, the agents may convert an 
appropriate amount of wealth proportional to their effi- 
ciency factor in their favour or against. Now, the model 
can be understood in terms of the general form of equa- 
tions: 



Wi{t + 1) = eiWi{t) -f ejWjit), 



Wj{t + !) = (!- e^)w^{t) + (1 - ej)wj{t), 



(23) 



where e^'s are quenched random numbers between and 
1 (randomly assigned to the agents at the beginning). 
Now the agents are supposed to make a choice to whom 
not to trade with. This option, in fact, is not unnatural 
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in the context of a real society where individual or group 
opinions are important. There has been a lot of works on 
the process and dynamics of opinion formations [J, |2l| 
in model social systems. In the present model it may be 
imagined that the 'choice' is simply guided by the rela- 
tive class index of the two agents. It is assumed that an 
interaction takes place when the ratio of two class fac- 
tors remain within certain upper limit. The requirement 
for interaction (trade) to happen is then 1 < ti/tj < r, 
where e^ > e^. Wealth distributions for various values of 
T are numerically investigated. Power laws in the tails of 
the distributions are obtained in all cases. In Fig.llUlthe 
distributions for t = 2 and r = 4 are shown. Power laws 
are clearly seen with an exponent, a = 3.5 (a straight line 
with slope around -3.5 is drawn) which means the Pareto 
index i/ is close to 2.5. It is not further investigated 
whether the exponent (a) actually differs in a significant 
way for different choices of r. It has been shown that 
preferential behaviour 20] generates power law in money 
distribution with some imposed conditions which allows 
the rich to get higher probability of getting richer. The 
rich is also favoured in a model with some kind of asym- 
metric exchange rules as proposed in 21] where a power 
law results in. The dice seems to be loaded in favour of 
the rich otherwise the rich can not be the rich! 
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FIG. 10: Distribution of individual wealth with selective in- 
teraction. Power law is evident in the log-log plot where a 
straight line is drawn with slope = -3.5 for comparison. 



VIII. MEASURE OF INEQUALITY 

Emergence of Pareto's law signifies the existence of in- 
equality in wealth in a population. Inequality or dis- 
parity in wealth or that of income is known to exist in 
almost all societies. To have a quantitative idea of in- 
equality one generally plots Lorenz curve and then cal- 
culates Gini coefficient. Here the entropy approach '39] is 
considered. The time evolution of an appropriate quan- 
tity is examined which may be regarded as a measure of 
wealth-inequality. 



Let us consider wi, W2, . . . , w^ be the wealths of N 
agents in a system. Let W = J2i=i ^i be the total wealth 
of all the agents. Now pi = Wi/W can be considered as 
the fraction of wealth the i-th agent shares. Thus each of 
Pi > and X]i=i P» = 1- Thus the set of pi, _P2, ■ ■ ■ , Pn 
may be regarded as a probability distribution. The well 
known Shannon entropy is defined as the following: 



N 



S = -y^pilnpi 



(24) 



From the maximization principle of entropy it can be 
easily shown that the entropy {S) is maximum when 



Pi =P2 



-Pn 



1 

iV' 



(25) 



giving the maximum value of S to be In N where it is 
a limit of equality (everyone possesses the same wealth) . 
A measure of inequality should be something which mea- 
sures a deviation from the above ideal situation. Thus 
one can have a measure of wealth-inequality to be 



Af 



N 



H ^\nN-S = \nN+^p,\nPi ^^p^\n{Np.i). (26) 
The greater the value of H, the greater the inequality is. 
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FIG. 11: Comparison of time evolution of the measure of 
inequality [H) in wealth for different models. Each 'time 
step' it) is equal to a single interaction between a pair of 
agents. Data is taken after every 10* time steps to avoid 
clumsiness and each data point is obtained by averaging over 
10^ configurations. Y-axis is shown in log-scale to have a fair 
comparison. 

It is seen that the wealth exchange algorithms are so 
designed that the resulting disparity or variance (or mea- 
sure of inequality), in effect, increases with time. Wher- 
ever power law in distribution results in, the distribu- 
tion naturally broadens which indicates that the variance 
(cr^) or the inequality measure [H in eqn. (|26|l ] should in- 
crease. In the Fig. ^2 and in Fig. E| time evolution of 
inequality measure H and variance u^ respectively are 
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than that is obtained in such a way. The resulting distri- 
bution of wealth is plotted in a log-log scale in Fig. 1131 



FIG. 12: Evolution of variance (a^) with 'time' (t) for dif- 
ferent models. F-axis is shown in log scale to accommodate 
three sets of data in a same graph. Data is taken after ev- 
ery 10^ time steps to avoid clumsiness and each data point is 
obtained by averaging over 10'^ configurations. 



plotted with time for three models to have a compar- 
ison. It is apparent that the measure of inequality in 
steady state attains different levels due to different mech- 
anisms of wealth exchange processes, giving rise to dif- 
ferent power law exponents. The growth of variance is 
seen to be different for different models considered, which 
is responsible for power laws with different exponents as 
discussed in the text. The power law exponents (a) ap- 
pear to be related to the magnitudes of variance that are 
attained in equilibrium in the finite systems. 



IX. DISTRIBUTION BY MAXIMIZING 
INEQUALITY 



It is known that probability distribution of wealth of 
majority is different from that of handful of minority 
(rich people). Disparity is more or less a reality in all 
economies. A wealth exchange process can be thought 
of within the present framework where the interactions 
among agents eventually lead to increasing variance. It is 
numerically examined [l9| whether the process of forcing 
the system to have ever increasing variance (measure of 
disparity) leads to a power law as it is known that power 
law is usually associated with infinite variance. Evolution 
of variance, a^ — (w^) — (w)^ is calculated after each in- 
teraction in the framework of pure gambling model [the 
pair of equations Q] and it is then forced to increase 
monotonically by comparing this to the previously calcu- 
lated value (the average value w is fixed by virtue of the 
model). This results in a very large variance under this 
imposed condition. The inequality factor H also like- 
wise increases monotonically and attains a high value. A 
power law distribution is obtained with the exponent, a 
close to 1. None of the available models does bring out 
such a low value of the exponent. The variance in any of 
the usual models generally settles at a level much lower 
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FIG. 13: Wealth distribution by maximizing the variance in 
the pure gambling model. Power law is clearly seen (in the 
log-log plot) and a straight line is drawn with slope — -1.1 to 
compare. 

Power law, however, could not be obtained by the same 
way in the case of a non-conserved model like the follow- 
ing: Wiit -|- 1) = 'Wi{t) ± S, where the increase or decrease 
(d) in wealth (w) of any agent is independent of any other. 

It has also been noted, considering some of the avail- 
able models, larger the variance, smaller the exponent 
one gets. For example, the variance is seen to attain 
higher values (with time) in the model of distributed 
(random) saving propensities [l^ compared to the model 
of selective interaction [l9j and the resulting power law 
exponent a is close to 2.0 in the former case whereas it 
is close to 3.5 in the later. In the present situation the 
variance attains even higher values and the exponent a 
seems to be close to 1, the lowest among all. 



X. CONFUSIONS AND CONCLUSIONS 

As it is seen, the exchange of wealth in various 
modes generates a wide variety of distributions within 
the framework of simple wealth exchange models as dis- 
cussed. In this review, some basic structures and ideas 
of interactions are looked at which seem to be funda- 
mental to bring out the desired distributions. In this 
kind of agent based models (for some general discussions, 
see '3.(fl) the division of labour, demand and supply and 
the human qualities (selfish act or altruism) and efforts 
(like investments, business) which are essential ingredi- 
ents in classical economics are not considered explicitly. 
What causes the exchange of wealth of a specific kind 
among agents is not important in this discussion. Mod- 
els are considered to be conserved (no inflow or outflow 
of money/ wealth in or from the system). It is not es- 
sential to look for inflation, taxation, debt, investment 
returns etc. of and in an economic system at the outset 
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for the kind of questions that are addressed here. The 
essence of complexities of interactions leading to distri- 
butions can be understood in terms of the simple (mi- 
croscopic) exchange rules much the same way the simple 
logistic equations that went on to construct 'roads to 
Chaos' and opened up a new horizon of thinking of a 
complex phenomenon like turbulence |37j . 

Some models of zero sum wealth exchange are exam- 
ined here in this review. One may start thinking in a 
fresh way how the distributions emerge out of the kind 
of algorithmic exchange processes that are involved. The 
exchange processes can be understood in a general way 
by looking at the structure of associated 2x2 transi- 
tion matrices. Wealths of individuals evolve to have a 
specific distribution in a stead state through the kind 
of interactions which are basically stochastic in nature. 
The distributions shift away from Boltzmann-Gibbs like 
exponential to Gamma type distributions and in some 
cases distributions emerge with power law tails known 
as Pareto's law {P{w) oc w^"). It is also seen that the 
wealth distributions seem to be influenced by personal 
choice. In a real society, people usually do not inter- 
act arbitrarily rather do so with purpose and thinking. 
Some kind of personal preference is always there which 
may be incorporated in some way or other. Power law 
with distinctly different exponent (a = 3.5, Pareto ex- 
ponent 1^ = 2.5) is achieved through a certain way of 
selective interaction. The value of Pareto index i^ does 
not correspond to what is generally obtained empirically. 
However, the motivation is not to attach much impor- 
tance to the numerical value at the outset rather than 
to focus on the fact of how power laws emerge with dis- 
tinctly different exponents governed by the simple rules 
of wealth exchange. 

The fat tailed distributions (power laws) are usually 
associated with large variance, which can be a measure 
of disparity. Economic disparity usually exists among a 
population. The detail mechanism leading to disparity 
is not always clear but it can be said to be associated 
with the emergence of power law tails in wealth distri- 
butions. Monotonically increasing variance (with time) 
can be associated with the emergence of power law in in- 
dividual wealth distributions. The mean and variance of 
a power law distribution can be analytically derived [3| 
to see that they are finite when the power law exponent 
a is greater than 3. For a < 3, the variance diverges 
but then the mean is finite. In case of the models dis- 
cussed here in this review, mean is kept fixed but large or 
enhanced variance is observed in different models when- 
ever there results in a power law. It remains a question 
of what can be the mechanisms (in the kind of discrete 
and conserved models) that generate large variance and 
power law tails. Large and increasing variance is also 
associated with lognormal distributions. A simple multi- 
plicative stochastic process like w{t+l) = e{t)w{t) can be 
used to explain the emergence of lognormal distribution 
and indefinite increase in variance. However, empirical 
evidence shows that the Pareto index and some other 



appropriate indices (Gibrat index, for example), gener- 
ally dwindle within some range jSS] indicating that the 
variance (or any other equivalent measure of inequality) 
does not increase forever. It seems to attain a saturation, 
given sufficient time. This is indeed the case the numer- 
ical results suggest. Normally there occurs simultaneous 
increase of variance and mean in statistical systems (in 
fact, the relationship between mean and variance goe s by 
a power law as a^ ex w know as Taylor's power law |39j 
as curiously observed in many natural systems). In this 
conserved model the mean is not allowed to vary as it is 
fixed by virtue of the model. It may be the case that a^ 
then has to have a saturation. The limit of cr^ is tested 
through an artificial situation where the association of 
power law with large variance is tested in a reverse way. 

Understanding the emergence of power law |7J, Ig it- 
self has been of great interest for decades. There is usu- 
ally no accepted framework which may explain the origin 
and wealth of varieties of its appearance. It is often ar- 
gued that the dynamics which generate power laws is 
dominated by multiplicative processes. It is true that 
in an economy wealth (or money) of an agent multiplies 
and that is coupled to the large number of interacting 
agents. The generic stochastic Lotka-Volterra systems 
like Wi(t+1) = ewi{t)+dw{t)—bwi{t)w{t) have been stud- 
ied p34l l4Clj| to achieve power law distributions in wealth. 
However, these kinds of models are not discussed in this 
review as the basic intention had been to understand the 
ideas behind models of conserved wealth which the above 
is not. 

In a twist of thinking, let us imagine a distribution 
curve which can be stretched in any direction as one 
wishes to have, keeping the area under this to be in- 
variant. If now the curve is pulled too high around the 
left then the right hand side is to fall off too quickly, ex- 
ponential decay is a possible option then. On the other 
hand, if the width of it is to be stretched too far (dis- 
tribution becomes fat) at the right hand side, it should 
then decay fairly slowly giving rise to a possible power 
law fall at the right end while keeping the area under the 
curve preserved. What makes such a stretching possible? 
This review has been an attempt to integrate some ideas 
regarding models of wealth distributions and to reinvent 
things with a fresh outlook. In the way, some confu- 
sions, conjectures and conclusions emerged where many 
questions possibly have been answered with further ques- 
tions and doubts. At the end of the day, the usefulness 
of this (review) may be measured by further curiosities 
and enhanced attention on the subject if at all this may 
generate. 
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